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Abstract. We present the first CCD photometric observations of the northern open cluster NGC 2126. Data were taken on 
eight nights in February and December 2002 with a total time span of ~57 hours. Almost 1000 individual V-band frames were 
examined to find short-period variable stars. We discovered six new variable stars, of which one is a promising candidate for an 
eclipsing binary with a pulsating component. Two stars were classified as 6 Scuti stars and one as Algol-type eclipsing binary. 
Two stars are slow variables with ambiguous classification. From absolute V(RI)c photometry we have estimated the main 
characteristics of the cluster: m - M = 1 l ra ± 0?5, E(V - 1) = 0M ± OF 1, E(V - R) = 0TO8 ± W06 (E(B - V) = 0T2 + 0™ 15) 
and d = 1.3 ± 0.6 kpc. Cluster membership is suggested for three variable stars from their positions on the colour-magnitude 
diagram. 
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1. Introduction 

Variable stars in clusters are crucial tools of stellar astro- 
physics: the application of stellar evolutionary theories via 
isochrone fitting of the colour-magnitude diagrams yields tem- 
perature, luminosity and age values for the member stars, 
which in turn draw strong constraints on pulsational properties 
(e.g. Frandsen & Arentoft 1998) or binary evolution (Rucinski 
et al. 1996). Early research of 6 Scuti pulsation in open clus- 
ters (Breger 1972) and its revival with the wide application 
of the CCD technique (for a review of the current state see 
Rodriguez & Breger 2001) have recently been extended toward 
the new class of y Dor variables (Handler 1999), which have 
been found in a number of young and intermediate-age open 
clusters (Arentoft et al. 2001, Kim et al. 2001, Choo et al. 2003 
and references therein). Clusters of different parameters can 
help map out the dependence of pulsational properties on age 
and metallicity, thus allowing better understanding of physical 
mechanisms driving the pulsation. On the other hand, detached 
eclipsing double-lined spectroscopic binaries can serve as very 
accurate distance indicators (e.g. Thompson et al. 2001), with 
the potential of feeding back to the isochrone fitting method it- 
self. This work aims to contribute to these issues with the first 
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CCD observations of the northern open cluster NGC 2126, with 
particular emphasis on its variable star content. 

NGC 2126 (= C0559+499, a 2 ooo=06 ,1 02'?55, 
<52000=+49°52', Triimpler class III 2 m :b) is a moder- 
ately rich typical galactic cluster with several dozens of 
members scattered in a region of 5-6' (Lynga 1987) in the 
constellation Auriga. The only previous photometric study 
was presented by Cuffey (1943), who made blue and red 
(A.A. 4300-6200) photographic observations. He estimated the 
distance to NGC 2126 as 950 pc, based on a comparison to 
M 35. The WEBDA catalog 1 reflects this lack of data, as no 
other source is listed for NGC 2126. This neglect turned our 
attention toward NGC 2126 and the present study presents 
an analysis of 8 nights of observations obtained in 2002. The 
measurements and data reduction are described in Sect. 2. 
Cluster parameters and the six new variable stars are discussed 
in Sect. 3, and a brief summary is given in Sect. 4. 



2. Observations and data reduction 

CCD V(RI)c observations were carried out on five consecu- 
tive nights in February 2002 and three nights in December 
2002. We used the 60/90/180 cm Schmidt telescope mounted 



1 http://obswww.unige.ch/webda 



2 



A. Caspar et al.: NGC 2126: physical parameters and variable stars 



Table 1. The journal of observations. 



Date 


V 


R c 


7 C 


length 
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frames 
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(hours) 


2002 Feb. 1/2 


83 


2 


3 


8.5 


2002 Feb. 2/3 


108 
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8.8 


2002 Feb. 3/4 


142 
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8.7 


2002 Feb. 4/5 


134 
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2002 Feb. 5/6 
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2002 Dec. 23/24 
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Total: 8 nights 
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Fig. 1. The observed CCD field (29' x 18') of the open cluster 
NGC 2126. The new variable stars are shown with open circles. 
Symbol sizes are proportional to the stellar brightnesses. The 
largest concentration of stars is about 8' in diameter, centred at 
(750,550). North is up, east is to the left. 



at the Piszkesteto Station of the Konkoly Observatory. The de- 
tector was a Photometries AT-200 CCD camera (1536x1024 
pixels, KAF-1600 chip with UV-coating). The image scale was 
l'.'l/pixel, giving a 29' x 18' field of view. The observations 
consisted of V-band time-series observations (mostly 180 s ex- 
posures) and calibrated V(RI)c photometric observations (from 
60 s to 600 s exposures). The full observing log is presented 
in Table ^ The observed field and the new variable stars are 
shown in Fig.^ 

The image reduction was done with standard tasks in 
IRAF 2 . For flat field corrections, we used sky flat images 
taken during the evening and morning twilights. We performed 
psf-photometry with the daophot package in IRAF using the 
Moffatt point-spread function, which was found to be well- 
suited for fitting the slightly distorted stellar profiles of the used 
instrument (Kiss et al. 2001). 

For photometric calibrations, we observed standard stars in 
the open cluster M 67 from the list of Chevalier & Ilovaisky 

2 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



(1991). The instrumental magnitudes were transformed with 
the photcal package, according to the following standard trans- 
formation equations: 

V = v-0.136X + 0.126(V-i?)- 5.180, <x = OTO (1) 
(V-I) = 1.028((v-0 -0.079Z) + 0.348, cr=0?02 (2) 
(V-R) = 0.827((v - r) - 0.034X) - 0.007, cr = 0™04 (3) 

where the symbols have their usual meanings. For the cal- 
ibration we used 19 stars of the 29 listed by Chevalier & 
Ilovaisky (1991). The colour ranges were -fT03...0? 1 64 and 
-0?07...1™21 for V - R and V-I, respectively, so that tran- 
formations of the reddest stars in the field (about 10% of 
the sample) were extrapolations. Consequently, there might be 
some systematic distorsion in the reddest part of the colour- 
magnitude diagram. 

The transformation of the cluster was done in two steps. 
First, we calibrated non-saturated stars on three images (one 
for each filter) confined by consecutive M67 observations on 
February 2, 2002. Then we co-added a number of V band im- 
ages and all (RI)c frames and measured fainter stars relative 
to previously calibrated local standards. We estimate the final 
photometric error as +0^05, which is likely to be an overesti- 
mate for the brightest stars. 

Time-series data were reduced in two separate subsets (5 
nights in February and 3 nights in December), because the 
follow-up observations in December, 2002 were obtained un- 
der a different instrumental setup (the CCD imager was rotated 
by 90°, so that a certain field of view was excluded from the 
measurements). Within the subsets, CCD frames were matched 
by the task imalign. Since the cluster is not too crowded, instru- 
mental magnitudes could be easily identified based on the X-Y 
pixel coordinates. Stars lying closer than 30 pixels to the image 
edge were rejected. Differential light curves were calculated 
with the ensemble normalization technique (Gilliland & Brown 
1988), for which dozens of stars of medium brightness were se- 
lected. We could identify almost 800 stars in the February run 
and that was the light curve dataset in which we searched for 
variable stars. The search was performed by checking the light 
curve statistics and by visual inspection of the suspect time- 
series data. We calculated the scatter of individual light curves 
for every star and plotted it against the apparent brightness. 
Outstanding points referred to either variable stars or constant 
stars with psf-fitting difficulties (e.g. crowding). In this way we 
found six stars to be intrinsic variables. In December, 2002, we 
focused on these six variables and the ensemble comparison 
stars, so that it is possible that there are unnoticed slow and/or 
long-period variables in the cluster field. 

The light curves were further analysed to search for 
possible periodicities. For this, we calculated Fourier spectra 
with Period98 (Sperl 1998). For the eclipsing binaries, the 
phase dispersion minimisation (PDM, Stellingwerf 1978) 
provided a better period searching method. Our original 
data are available electronically at the CDS via anonymous 
ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.£r/cgi-bin/qcat?J/A+A/. . ./ 
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Fig. 2. The colour-magnitude diagrams of NGC 2126. Top row : 
full field of view; bottom row: the central field of 8' diameter. 
The solid lines show the log t — 9.1 isochrone. 



3. Results 

3.1. The colour-magnitude diagram 

The physical parameters of the cluster were estimated by fitting 
isochrones to the colour-magnitude diagram (CMD). In order 
to separate the cluster's stars from other field stars, we have 
examined proper motions taken from the USNO B1.0 catalog 
(Monet et al. 2003). We found that of the 800 stars, almost 
200 have proper motions detected in the USNO plate material. 
However, those stars are evenly distributed over the field with 
no concentration around the cluster, so that NGC 2126 turned 
out to be in the background. Consequently, we removed all stars 
from the CMD with non-zero proper motions and the resulting 
CMDs are shown in the top row of Fig. [2 This procedure has 
been checked by two nearby galactic fields, offset by about 1 
degree north and south. A comparison of proper motion his- 
tograms showed that we have indeed selected the galactic fore- 
ground. 

In order to decrease background contamination, we kept 
only the inner 8' of the cluster for the isochrone fitting (bot- 
tom row in Fig. 0. This region contains 103 stars between 
V = 1 1?3 - 18™5, of which the majority seems to be easily dis- 
tinguishable from the background and some weak foreground. 
We assumed solar composition and included reddening deter- 



mination in the fitting procedure. The isochrones were taken 
from Bertelli et al. (1994) and shifted individually to match 
the main sequence, turn-off point and red giant positions. We 
found that the overall shape of the CMD is well reproduced 
with isochrones of log t ranging from 9.0 to 9.3, E(V - I) be- 
tween 0?53 and (T33, E(V-R) between 07> 14 and (F03 and dis- 
tance moduli between 11™5 and 107M-. We plot the logf = 9.1 
isochrone in Fig. [2] showing the "best" fit, but the difference for 
the other isochrones is almost negligible. Changing the com- 
position makes the parameter ranges even wider, thus the es- 
timated parameter errors are quite large and a further spectro- 
scopic study addressing spectral types and chemical composi- 
tion is desirable. 

In this paper we adopt the following physical parameters: 
m-M= l¥?0±O?6,E(y-I) = 07 1 4 + 0?1,£(V-/?) = 07K)8± 
0^06. The reddenings were converted to E(B - V) using the 
coefficients listed in Rieke & Lebofsky (1985): E(V-R)/E(B- 
V) = 0.78, E(V - I)/E(B - V) = 1.60. The result is E(B - V) = 
0^2 + 07 1 15, and the distance to the cluster (assuming Ay « 
3.1 xE(B-V)) is 1.3 + 0.6kpc. 

3.2. New variable stars 

Nightly light curves of the new variable stars are shown in 
Fig. |3] (numbered in order of increasing right ascension). 
Differential light curve data were shifted so that the mean val- 
ues match the apparent magnitudes of the stars in the CMD. 
Therefore, the absolute V magnitudes in Fig. [3] are somewhat 
uncertain, depending on the magnitude range of the variables. 

The light curve shapes suggest the following classifica- 
tions: VI and V2 are slow variables of ambiguous nature; V3 
and V5 are short-period variables, most likely of 6 Scuti type; 
V4 is an Algol-type eclipsing binary with only one observed 
minimum, while the light curve of V6 is a mixture of 5 Scuti- 
like oscillations and Algol-like eclipses. 

3.2.1 . V1 and V2: slow variations 

Although these two stars showed clear variability, our dataset 
is too short to determine reliable periods. For VI, both PDM 
and Fourier analysis suggested P= 1.64470 d or half that value, 
0.82235 d (the uncertainties are a few in the last digit). We plot 
the resulting phase diagrams in Fig.|4](the epoch of maximum 
is HJD 2452307.48). 

With no colour light curve, the classification is uncertain. 
For 0.822235 d, one can speculate on possible y Dor-like oscil- 
lations. These stars are early F-type main-sequence or subgiant 
stars with periods 0.4-3.0 d and V-band amplitudes less than 
Wl, pulsating in non-radial gravity modes. According to the 
definition by Handler (1999), the y Dor instability strip cov- 
ers the 7200-7700 K range on the zero-age main sequence 
(ZAMS) and the 6900-7500 K range one magnitude above the 
ZAMS. Assuming cluster membership, the dereddened colours 
of VI are (V-R) = 0^24 and (V-/) = 0>M8. Standard tabula- 
tions give an effective temperature of 6500-7000 K (Cox 2000), 
which marginally supports the classification as a y Dor star. 
Also, the cluster membership would imply an absolute mag- 
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Fig. 3. Individual light curves of the six new variable stars. The solid lines denote Fourier fits for V3 and V5. 



nitude My=3™0; the temperature and absolute magnitudes are 
consistent with the coolest y Dor variables. In this case, the 
apparent amplitude modulation is due to the presence of other 
excited modes of pulsation and a proper asteroseismological 
analysis would require a much longer dataset. 

On the other hand, for P= 1 .6447 d, the double-peaked light 
curve is similar to those observed for spotted active binaries 
of RS CVn subtype. For instance, recent data collection of 
Cutispoto et al. (2003) contains nice examples of very simi- 
lar double-peaked curves of well-known active stars (BI Cet, 
BC Phe). Based on the present dataset, it is impossible to dis- 
tinguish between the two classifications. The fact that VI lies 
quite far from the cluster centre (about 10') weakens the possi- 
bility of membership and consequently makes the temperature 
and absolute magnitude estimates less reliable. 

For V2, the data can be equally well folded with periods 
around 0.5 d and 1 d, but neither period gives a continuous 
phase diagram with no gap. Therefore, we cannot decide the 
cause of the low-amplitude variability of this star (the total 
range is about CT08). 

3.2.2. V4: an Algol-type binary 

We observed only one minimum of the star (at HJD 
2452308.387), so that no firm conclusion can be drawn on its 



period. However, if we plot only the February subset in one 
light curve (Fig. [5}, the shape suggests a clearly visible reflec- 
tion effect (about Assuming symmetry around the un- 
observed secondary minimum, a tentative period of ~3 days 
can be estimated (i.e. the next primary minimum happened 
immediately before or after the last night of observations in 
February). 

3.2.3. V3 and V5: multiply periodic oscillations 

Both stars showed rapid oscillations with full amplitudes of a 
few tens of millimags, characteristic of 6 Scuti variables. We 
have determined the pulsational frequencies with Fourier anal- 
ysis consisting of iterative prewhitening steps (see Fig. [6). 

We could identify two frequencies in both light curves (see 
Table EJ with S/N ratios (Breger et al. 1993) larger than 4. 
For V5, the low-frequency component refers to slight night-to- 
night variations of the mean brightness level, which we suspect 
were caused by a closely separated field star (see Fig.0 lying 
~ 7" from V5. It is fainter by about 4 mag, so that it may affect 
profile fitting to some small extent. Nightly variations of the 
seeing may also introduce variable contribution from this star. 

In summary, V3 and V5 are probably multiply periodic pul- 
sating 5 Scuti variables. The frequency ratios (V3: fi/fi = 
0.81, V5: fi/fi = 0.94) suggest non-radial modes of pulsa- 
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Fig. 4. Phase diagrams for VI with two equally possible peri- 
ods. 
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Fig. 5. The 2002 February data for V4. 



tion for both stars. Small deviations from the light curve fits 
(Fig. |3 and residual structures of the prewhitened spectra (es- 
pecially for V3) suggest the probable existence of more fre- 
quencies. Assuming their cluster membership implies the fol- 
lowing dereddened colours and absolute magnitudes: V3 — 
(V - R) Q = 0^25, (V - 7)o = 0?5, M v « 4 m 4; V5 — 
(V - K) = 0^23, (V - 7)o = 0M3, My * 2™0. The latter values 
are consistent with typical parameters of 6 Set stars (Rodriguez 
& Breger 2001), we therefore conclude V5 is likely to be a 
member of the cluster. Moreover, its proximity to the cluster 
centre supports this conclusion. On the other hand, V3 seems to 
be a background object, since its absolute magnitude, assuming 
cluster membership, places it well below the main sequence. 

3.2.4. V6: an eclipsing binary with a pulsating 
component? 

The most intriguing star is V6, which shows complex light vari- 
ations. We observed well-defined minima - three deeper and 
two shallower - characteristic of an eclipsing binary with com- 
ponents of different temperatures. Moreover, the star showed 
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Fig. 6. Fourier analysis of V3 and V5. 



Table 2. The results of the period analyses. 
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12.14 
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steady oscillations outside eclipses, with amplitude and cycle 
length characteristic of 5 Scuti-like pulsation (about 0^05 and 
0.13 d). The symmetric stellar profile, even on the best CCD 
images, excluded the possibility of two unresolved stars, so that 
we conclude V6 is an eclipsing binary with at least one pulsat- 
ing component. Such systems are spectacular targets for astero- 
seismology, because independent determination of the physical 
parameters (mass, radius, temperature), applying binary star as- 
trophysics, gives strong constraints on the possible mode iden- 
tification (e.g. Mkrtichian et al. 2002, Kim et al. 2002a). In ad- 
dition, eclipses of a pulsating star give spatial resolution across 
the stellar surface and thereby provide the possibility for mode 
identification (e.g. Reed et al. 2002). 

The period analysis resulted in P or b = 1.17320(3) d. The 
data were folded with this period and the phase diagram is 
shown in Fig.|7] The depth is (F2 and 0? 12 for the primary and 
secondary minima, respectively. The oscillations are clearly 
visible outside minima and we note that the scatter of the phase 
diagram after folding only the February data was fairly small. 
This indicated that the oscillations were quite coherent within 
5 consecutive nights in February. Their amplitude was about 
0"05 associated with a period P pu i ^0.13 d. 

We also performed a separate period analysis of data from 
which eclipse minima were excluded. This showed the oscil- 
lations to be surprisingly stable. They have a period f pu i = 
0.12936 d and the phase diagram shows remarkable small 
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Fig. 7. Phase diagram for V6. 
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Fig. 8. Outside-eclipse data folded with P pu i = 0.12936 d. The 
zero phase corresponds to a primary eclipse. 
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Fig. 9. A low -resolution spectrum obtained on April 10, 2003. 



nearby objects residing south-east to it. In this way, we mini- 
mized the possibility of light from the nearby objects contam- 
inating the V6 spectrum. The spectrum shows prominent hy- 
drogen absorption lines accompanied by the sodium D-doublet 
(Fig. [5}. The spectral lines and the continuum shape suggest 
an F-type star, typical for a S Scuti star, so that the oscillations 
may be attributed to 6 Scuti-like pulsations. 

Out of the six variables, only V6 has non-zero proper mo- 
tion (fi a cos <5 = —14 mas/year, fig = 10 mas/year), thus cluster 
membership can be excluded. 

In the literature, only a few similar system are reported, 
making V6 a very interesting variable star. Rodriguez & Breger 
(2001) listed nine 6 Scuti variables in eclipsing binaries, of 
which RZ Cas (Ohshima et al. 2001) and AB Cas (Rodriguez 
et al. 1998) are the best studied. Compared to the known ex- 
amples, V6 is different in two important respects: i) the relative 
amplitude ratio (AV pr i m /AV pu i ~ 4) is much smaller than that 
of any other similar stars (the next is V577 Oph, with a ra- 
tio of 13); ii) the period ratio {P 01 b/P P ui ~ 9) is much smaller 
than usual - from Table 3 in Rodriguez & Breger (2001), only 
WX Eri has a ratio smaller than 20 (5.005). 

Very recently, several new examples of pulsating compo- 
nents in eclipsing binaries have been found (RX Hya, Kim et 
al. 2002b, 2003; AB Per, Kim et al. 2002c, 2003; TW Dra, 
Kim et al. 2003; AS Eri, Gamarova et al. 2000), but all are typ- 
ical low-amplitude rapid pulsators not resembling V6. In that 
sense, V6 is a unique candidate for an eclipsing binary with a 
pulsating component. Future spectrophotometric observations, 
accompanied by radial velocity measurements, are needed to 
settle a number of questions regarding V6. In particular, are the 
cyclic changes caused by stellar pulsation? Which component 
is pulsating? What is the true amplitude of the oscillations, af- 
ter removing the second light from the light curve (could this 
be the first high-amplitude 6 Scuti star in an eclipsing binary)? 
Is there indeed resonance between orbital motion and pulsa- 
tion? If so, are the oscillations induced by tidal forces (Willem 
& Aerts 2002)? 

4. Conclusions 

This paper presents the first CCD photometric study of the 
northern open cluster NGC 2126. The V(RI)c observations re- 
vealed the main characteristics of the cluster and led to the dis- 
covery of six new variable stars. The locations of the variables 
on the CMD (Fig.llQI suggest membership for VI, V4 and V5; 
V2 and V3 are most likely non-member stars, while V6 is a 
foreground object. The most important results of this research 
are the following: 



scatter (Fig. The subsequent prewhitening steps indicated 
a low-frequency component, approximately 2/T bs, and fi a 
fpui - /orb- Another interesting result is that /p U i//oi-b = 9.07, 
suggesting that there might be a 1 :9 resonance between the or- 
bital motion and pulsation. 

We have also observed V6 with the 40" telescope at the 
Tel- Aviv University Wise Observatory on 2003 April 10. The 
spectrograph 10" long-slit was aligned at PA=133.8° in order 
to maximize the separation between the V6 spectrum and the 



1. We have estimated important physical parameters of the 
cluster with standard photometric methods for the first 
time. 

2. We have discovered six variable stars in a field of 29' x 18' 
centred on the cluster, which is about 8' in diameter. Two 
are pulsating stars, one is an eclipsing binary with unknown 
period, one showed a combined light variation of eclipses 
and pulsations, while two stars have ambiguous classifica- 
tion. 
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Table 3. The basic data of new variable stars. Except V6, coordinates and identifications were taken from the Guide Star Catalog, 
Version 2.2 (STSci 2001). 
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Fig. 10. Variable stars on the colour-magnitude diagram ex- 
cluding foreground stars (except V6). 



3. The most interesting variable star is V6, which is a promis- 
ing candidate for an eclipsing binary with a pulsating com- 
ponent. The period ratio Porb/^pui = 9.07 may indicate a 
resonance between the orbital motion and pulsation. 

We summarize the basic data of the variable stars in Table [3] 
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